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SECT  I ON  1 
INTRODUCTION 

This  paper  is  basically  an  elaboration  and  an.  extension 
of  some  of  our  earlier  wor*k  that  has  already  been  published  [ 1 ] ; 
however,  there  is  an  important  difference  between  that  work 
and  this.  Aside  from  the  extension,  the  difference  is  that  at 
this  t,ime  we  are  fortunate  to  have  high  quality  laboratory 
data  available  [2]  for  comparison  with  our  calculations.  We 
intend  to  take  advantage  of  the  degree  of  agreement  of  our 
calculations  with  these  data  to  establish  confidence  in  the 
credibility  of  our  calculations  for  situations  where  no 
laboratory  data  are  available. 

Our  claim  is  that  an  underlying  assumption  often  used  to 
justify  the  gross  modeling  of  the  physical  features  of  metallic 
enclosures  is  in  error.  In  order  to  clarify  the  relevance  of 
the  claim,  the  term  "metallic  enclosures"  is  meant  to  include 
structures  such  as  aircraft,  missiles,  ships,  and  tanks.  The 
specific  assumption  challenged  is  that  if  the  wavelength  of 
an  incident  monochromatic  plane  wave  is  long  compared  to  the 
dimensions  of  a metallic  structure,  then  the  detailed  features 
of  the  structure  are  not  important  for  the  subsequent  modeling 
of  the  structure  for  EMP-related  experiments  or  calculations. 

It  is  conjectured  that  the  historical  basis  for  this  assump- 
tion is  its  validity  for  far  zone  calculations  or  measurements. 
The  argument  relating  this  frequency  domain  assumption  to  time 
domain  EMP  considerations  is  that  the  wavelengths  in  most  of 
a representative  EMP  spectrum  are  long  compared  to  many  physical 
features  of  interest  and  we  are  not  taking  issue  with  this 
observation. 

Our  evidence  to  support  our  claim  is  based  on  external 
coupling  analyses  and  measurements,  and  in  presenting  the 


1 


evidence  we  will  focus  our  attention  on  the  determination  of 
surface  current  densities.  The  reason  for  focusing  on 
densities  is  that  it  is  these  quantities  that  are  a necessary 
part  of  the  external  coupling  information  required  by  commonly 
employed  deliberate  antenna  and  aperture  analyses  for  the 
subsequent  calculation  of  voltages  and  currents  driving 
mission-critical  subsystems  contained  within  the  enclosure. 
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SECTION  2 

OUTLINE  OF  EVIDENCE  AND  ARGUMENTS 


We  intend  to  establish  confidence  that  our  finite  circular 
cylinder  MF IE  code  can  give  reliable  information  where  no 
measurements  are  readily  available  to  us.  In  order  to  do  this, 
we 

•'  compare  our  current  density  calculations  with  those 
measured  in  a laboratory  experiment 

• argue  that  our  code  should  perform  at  least  as  well 
at  lower  frequencies  corresponding  to  the  bulk  of  the 
energy  for  a typical  spectrum  of  an  EMP . 

We  will  then  present  current  density  results  obtained  by 
running  our  finite  circular  cylinder  code  for  a range  of  low 
frequencies  and  show  that  an  appropriate  magnetostatic  current 
density  solution  describes  the  dominant  behavior  of  the  current 
density  for  an  extended  range  of  frequencies.  Next  we  will 
explain  how  this  magnetostatic  effect  invalidates  commonly 
employed  external  coupling  approaches  that  assume  a uniform 
circumferential  distribution  for  the  current  density  in  order 
to  determine  this  density  simply  by  dividing  a calculated  bulk 
current  by  the  circumference.  In  regard  to  this  common 
assumption  we  will 

• show  that  the  density  calculation  becomes  increasingly 
more  inaccurate  as  the  commonly  stated  validity  condi- 
tions are  increasingly  better  satisfied 

• explain  how  the  physical  interpretation  of  a nonphysical 
situation  is  a likely  candidate  for  the  cause  of  the 
invalid  assumption 

• describe  the  impact  on  aircraft  stick  models  and 
certain  missile  models. 
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The  remaining  material  that  will  be  presented  is  new  in 
that  it  is  what  we  referred  to  as  the  extension  of  Reference  1. 
We  will  present  low  frequency  calculations  using  the  MFIE 
approach  for  the  current  density  induced  on  a finite  cylinder 
having  an  elliptic  cross-section.  We  then  compare  these 
finite  elliptic  cylinder  results  with  an  appropriate  magneto- 
static solution  and  again  demonstrate  that  the  dominant 
behavior  is  described  by  the  magnetostatic  solution  for  an 
extended  range  of  frequencies.  Using  the  magnetostatic  solu- 
tion, we  obtain  an  analytic  expression  for  the  error  that 
persists  at  low  frequencies  which  is  made  by  modeling  the 
cross-section  of  an  elliptic  cylinder  as  a circle.  The  fact 
that  there  is  any  long  wavelength  error  at  all  demonstrates 
that  long  wavelength  arguments  do  not  provide  justification 
for  gross  modeling. 

Finally,  we  present  results  based  on  symmetry  arguments 
to  qualitatively  demonstrate  the  distribution  of  current 
density  on  a complex  model  of  an  aircraft  at  low  frequencies. 
These  results  are  equivalent  to  a magnetostatic  solution,  and 
consequently  they  differ  from  what  would  be  obtained  by  a 
modeling  that  is  insensitive  to  magnetostatic  considerations. 
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SECTION  3 

COMPARISON  WITH  EXPERIMENTAL  DATA 

The  material  presented  in  the  first  five  figures  in  this 
paper  contains  material  that  was  presented  in  two  recent  reports 
[2,3].  The  curves  in  Figure  2 and  the  top  portions  of  Figures 
3,  4,  and  5 are  curves  that  were  traced  from  Xerox  copies  of 
curves,  presented  in  those  reports.  The  description  of  the 
experiment  pertaining  to  the  data  presented  in  these  figures 
is  also  redrawn  based  on  a figure  presented  in  Reference  2. 

The  intent  of  this  experiment  is  to  simulate  a monochromatic 
plane  wave  incident  on  a tube  having  a total  length  of  2h.  By 
referring  to  Figure  1,  we  can  see  that  the  angle  0 is  defined  so 
that  0°  corresponds  to  the  deep  shadow  region  and  180°  corresponds 
to  direct  illumination.  The  0 in  Figures  3,  4,  and  5 corresponds 
to  this  definition  and  z is  the  axial  distance  ranging  from 
0 at  the  ground  plane  to  h at  the  top  of  the  cylinder.  The 
quantities  1x^,1  , | Kri  | , and  Pz  plotted  in  these  figures  are  the 
magnitude  of  the  axial  component  of  the  current  density, 
magnitude  of  the  transverse  component  of  the  current  density, 
and  the  phase  of  the  axial  component  of  the  current  density. 

The  same  quantities  calculated  by  our  MFIE  computer  code  for 
the  current  density  induced  by  the  same  source  on  a flatly 
capped  cylinder  having  the  same  length  and  diameter  as  the 
tube  of  the  experiment  are  presented  in  the  bottom  half  of 
Figures  3,  4,  and  5.  The  scales  of  our  calculations  were 
adjusted  to  the  scales  of  the  experimental  data  by  using 
three  numbers,  a multiplicative  factor  for  each  set  of  the 
magnitude  comparisons  and  an  additive  factor  for  the  set  of 
phase  comparisons.  These  three  numbers  were  determined  by 
forcing  one  point  of  the  experimental  data  to  match  one  point 
of  our  calculated  data  on  only  one  curve  of  each  of  the  three 
sets  of  curves.  The  reason  we  include  Figure  2 in  this  paper 
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Figure  2.  Measured  Magnitude  and  Phase  of  Surface  Density  of  Outside  Axial 
Current  on  Tubular  Cylinder  with  Open  End,  Flat  and  Hemispherical  End  Caps. 
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Figure  3a.  Measured  Amplitude  of  Axial  Surface  Density  of 
Outside  Current  on  Tubular  Cylinder,  E-Polarization  (Large 
Outdoor  Ground  Screen) 


Figure  3b.  Calculated  Amplitude  of  Axial  Surface  Current 
Density  on  a Finite  Circular  Cylinder 
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ure  4a.  Measured  Amplitude  of  Surface  Density  of  Transverse  Current  on 
Tubular  Cylinder,  E-Polarization;  Normal  Incidence 
(Large  Outdoor  Ground  Screen) 
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Figure  4b.  Calculated  Amplitude  of  Transverse  Surface  Current  Density  on 

a Finite  Circular  Cylinder 
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Figure  5a.  Measured  Phase  of  Axial  Surface  Density  of  Outside  Current  on 
Tubular  Cylinder.  Normalization:  -0obs  (Large  Outdoor  Ground  Screen) 


Figure  5b.  Calculated  Phase  cf  Axial  Surface  Current  Density  on  a Finite 

Circular  Cyl inder 
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is  to  show  that  there  is  only  a minimal  measured  effect  of 
capping  the  tube,  thus  justifying  our  comparison  of  the  capped 
tube  calculations  to  the  uncapped  measurements.  It  should  be 
noted  that  our  comparisons  with  the  data  were  for  h = 36  cm 
while  the  data  in  Figure  2 correspond  to  h - 84  cm. 

Now  it  is  necessary  for  us  to  discuss  the  frequency  at 
which  we  made  the  comparison  between  our  calculations  and  the 
experimental  data.  The  normalized  value  kh  = 1.5  u determines 
the  frequency.  First  we  note  without  scaling  h to  missile  or 
aircraft  size  dimension  that  the  comparison  was  made  well 
beyond  the  primary  resonance  of  the  cylinder.  Next  we  mention 
that  if  h is  taken  to  be  in  the  10  to  20  m range,  the  frequency 
scales  to  the  20  to  10  MHz  range.  The  main  claim  of  this 
paper  deals  with  low  frequency  behavior  clearly  below  the 
frequency  at  which  we  made  our  comparison.  In  this  regard, 
we  claim  that  our  MFIE  approach  can  only  perform  better 
[ as  the  frequency  is  decreased.  Without  too  much  elaboration, 

the  basic  reason  is  that  two  types  of  zoning  requirements  had  to 
be  satisfied  to  obtain  the  agreement  presented  at  kh  = 1.5  it. 

One  requirement  is  frequency  independent  and  depends  only  on 
the  geometry  of  the  structure  while  the  other  is  the  usual 
frequency  dependent  requirement  that  no  dimension  of  a zone 
should  be  greater  than  a tenth  of  a wavelength.  The  agreement 
at  kh  = 1.5  it  insures  that  we  have  adequately  determined  the 
geometry  requirement  on  the  zoning  and  the  frequency  require- 
\ ment  becomes  less  restrictive  as  the  frequency  is  decreased. 
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SECTION  4 

LOW  FREQUENCY  FINITE  CIRCULAR  CYLINDER  DATA 

In  Figure  6 we  present  the  values  of  the  axial  component 
of  current  density  K^Och,  z,  0)  obtained  by  running  our  code 
for  the  relatively  low  frequencies  corresponding  to  kh  ranging 
from  0.025  to  0.5.  The  real  frequencies  for  these  values  of  kh 
can  readily  be  computed  from  the  approximate  formula 

f ^ 4.8  (kh/d)  MHz 

where  d is  the  value  of  h expressed  in  tens  of  meters.  The 
reason  for  expressing  f in  this  form  is  that  d will  typically 
range  between  1 and  2 since  h will  range  between  10  and  20  m 
for  typical  aircraft  and  missile  dimensions.  The  skin  depth 
corresponding  to  this  range  of  frequencies  can  readily  be 
found  in  graphical  form  for  different  metals  by  referring  to  Ref- 
erence 4.  In  addition  to  the  skin  depth  being  small  for  the  band  of 
frequencies  corresponding  to  the  data  presented  in  Figure  6, 
nonplanar  shielding  studies  as  discussed  in  Reference  1 show 
that  the  frequency  for  which  real  metal  can  be  considered 
perfectly  conducting  is  even  lower  than  what  is  indicated 
by  having  the  skin  depth  being  small  compared  to  wall 
thickness.  In  summary,  real  metal  used  for  aeronautical 
systems  behaves  as  though  it  were  perfectly  conducting  for 
the  data  presented  in  Figure  6. 

The  previous  discussion  regarding  whether  real  metal 
can  be  considered  perfectly  conducting  at  these  frequencies 
was  presented  in  order  to  make  the  applicability  of  magneto- 
statics more  palatable  since  it  assumes  a perfectly  conducting 
boundary  condition.  Displayed  in  Figure  6 is  an  appropriate 
magnetostatic  solution,  the  derivation  of  which  appears  in 
a later  section  of  this  paper.  It  is  important  to  note  that 
the  magnetostatic  solution  describes  either  the  dominant 
behavior  or  a significant  part  of  the  behavior  for  all  the 
data  displayed. 
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SECTION  5 

CONSEQUENCES  OF  LOW  FREQUENCY  FINITE  CIRCULAR  CYLINDER  DATA 

One  of  the  consequences  is  readily  discussed  by  notiny  the 
column  of  data  headed  by  K (old).  These  entries  are  values 
of  the  current  density  calculated  by  first  inteyratiny  our 
current  density  with  respect  to  if,  cor  respond  iny  to  inteyration 
with  respect  to  circumferential  arc  length,  in  order  to  obtain 
a bulk  current-  Next  we  divide  this  bulk  current  by  the  cir- 
cumference of  the  cylinder  and  denote  the  result  as  K (old). 

The  reason  for  obtaininy  this  quantity  is  that  it  corresponds 
to  a common  approach  for  obtaininy  this  quantity.  Specifically, 
it  is  often  bulk  current  that  is  directly  calculated  and  then 
current  densities  are  assumed  to  be  simply  derivable  by 
dividiny  the  bulk  current  by  a local  circumference.  It  is 
assumed  that  this  method  of  obtaininy  the  current  density 
becomes  an  increasingly  better  approximation  as  the  wavelenyth 
increases  relative  to  cross-sect  ion  dimensions.  By  comparing 
Kz  (old)  to  the  actual  values  contained  in  the  table,  we  see 
that  it  becomes  an  increasingly  worse  approximation  as  the 
wavelength  is  increased,  i.e.,  kh  decreases  for  a fixed  h/a 
ratio.  The  fact  that  the  approximation  qets  worse  as  the 
stated  conditions  for  its  validity  become  increasingly  better 
satisfied  indicates  a fundamental  error  in  interpreting  the 
physics  of  the  situation.  Before  completing  this  aspect  of 
our  presentation,  we  would  like  to  clarify  the  fact  that  we 
are  not  claiming  that  the  bulk  current  for  fat  cylinders  has 
been  incorrectly  calculated.  In  fact,  we  obtain  excellent 
agreement  with  the  results  presented  in  Reference  5 by  inte- 
grating our  current  density.  The  claim  is  that  a perfect 
knowledge  of  the  bulk  current,  no  matter  how  it  is  determined, 
cannot  be  used  to  obtain  current  density  by  dividiny  this 
quantity  by  the  circumference. 
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The  consequences  of  this  erroneous  assumption  is  that 
considerable  effort  has  been  devoted  toward  obtaining  bulk 
current  both  computationally  and  experimentally  when  it  is 
actually  current  density  that  is  required.  Until  recently, 
the  entire  intent  of  aircraft  stick  models  has  been  to  obtain 
bulk  current  and  then  infer  current  density  by  dividing  by 
the  circumference.  Regarding  certain  missile  models,  we  were 
involv.ed  in  calculating  the  bulk  current  induced  on  the 
various  stages  of  Minuteman  in  flight.  We  then  employed  the 
assumption  that  we  have  demonstrated  to  be  erroneous  and 
divided  the  bulk  current  by  the  local  circumference  in  the 
attempt  to  obtain  the  current  density  at  various  points  of 
entry . 

Returning  to  aircraft  stick  models,  a recent  modification 
has  been  to  add  a magnetostatic  solution  to  the  current  density 
calculated  in  the  old  manner  [6).  This  has  been  shown  to 
improve  the  calculation  in  that  the  results  are  in  better 
agreement  with  laboratory  measurements  after  the  addition. 

The  full  implications  of  this  modified  approach  deserve  to 
be  investigated  and  along  these  lines  several  issues  are 
brought  to  mind  by  our  experience.  One  issue  is  the  deter- 
mination of  the  conditions  under  which  the  total  current 
density  separates  into  the  two  separately  calculable  parts 
and  another  is  how  does  one  simply  obtain  a magnetostatic 
solution  for  a complex  aircraft  model. 


15 


SECTION  6 


POSSIBLE  EXPLANATION  OE  ORIGIN  OF  ERRONEOUS 
BULK  CURRENT  ASSUMPTION 


In  this  section  we  will  show  how  the  physical  interpre- 
tation of  a nonphysical  situation  is  a possible  cause  for  the 
commonly  employed  erroneous  assumptions.  We  will  start  with 
the  standard  solution  for  the  current  density  induced  on  an 
infinitely  long  perfectly  conducting  cylinder;  however,  our 
interpretation  of  this  equation  is  not  the  standard 
interpretation.  For  the  geometry  depicted  below 


for 

E.  = eik°*  2 
— inc 

the  current  density  induced  on  the  cylinder  Kz  is  given  by 


(ttZo/2)Kz  (u,0)  = u 


-1 


w , . . n 
(£  ^ 


n=0 


( i ) a n cos  n6 
(a) 


in,.,  )iFI“'e) 

where  a = k^a  = tt  , and  a^=l,  a^  = 2 n _>  1 . 


o a o n 

We  now  decompose  F(u,0)  into  a part  that  is  independent 
of  0 and  a remaining  part  that  does  depend  on  0.  That  is 


F ( a , 0 ) = Find(a)  + R(a,6) 


where  Fincj(“)  = (aH^^(a))  1 , 
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We  now  note  that  F^nc|(a)  approaches  infinity  as  a approaches 
zero  while  R(a,G)  is  bounded  for  all  a.  The  limit  corresponding 
to  a becoming  increasingly  small  is  the  limit  used  to  justify 
the  approximation  that  the  current  density  is  uniformly  distri- 
buted around  the  circumference  of  the  cylinder,  i.e.,  indepen- 
dent of  0.  The  following  inequality  based  on  the  described 

behavior  of  F.  ,(u)  and  R(u,0)  is  taken  as  the  mathematical 
1 nd 

statement  of  the  uniformity  assumption 

FindU)  >:>  R(a'0) 

as  a approaches  zero.  Our  point  is  that  this  is  a misleading 
and  physically  meaningless  inequality  when  interpreted  for 
cylinders  of  finite  length.  As  a general  rule  it  is  dangerous 
to  infer  physical  information  concerning  the  relative  largeness 
of  a quantity  going  to  infinity  without  investigating  the 
reason  for  this  quantity  approaching  infinity.  In  the  situa- 
tion presented  here,  by  starting  with  the  infinite  cylinder 
solution  it  was  never  possible  to  determine  the  effect  of 
relating  the  wavelength  to  the  length  of  the  cylinder.  To 
provide  a final  argument  to  show  that  the  relative  largeness 
of  should  not  be  used  to  infer  finite  length  cylinder 

behavior,  we  note  that  the  bulk  current  induced  on  the  infinite 
cylinder  is 

/2it 

Kz(a,0)  dO  = (2a/TTZQ)  Find  (u) 
o 

From  this  equation  we  see  that  the  fact  that  F^nd(a)  approaches 
infinity  for  small  a causes  1(a)  to  approach  infinity  for  small 
a,  while  at  the  same  time  it  is  well  known  that  1(a)  for  a 
cylinder  of  finite  length  approaches  zero  as  a approaches 
zero . 
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SECTION  7 

LOW  FREQUENCY  FINITE  ELLIPTIC  CYLINDER  DATA 

Figure  7 contains  the  values  of  the  axial  component  of 
the  induced  current  density  Kz(kh,  h/12,  <p ) obtained  by  running 
our  code  for  the  relatively  low  frequencies  corresponding  to 
kh  ranging  from  0.02  to  0.2.  The  real  frequencies  for  these 
values  of  kh  can  be  readily  determined  as  described  in  Section  4, 
and  the  discussion  concerning  the  applicability  of  perfectly 
conducting  boundary  conditions  for  real  metal  at  these 
frequencies  as  presented  in  that  section  is  also  relevant 
for  the  data  described  in  this  section.  Again,  it  is  important 
to  note  that  the  magnetostatic  solution  describes  either  the 
dominant  behavior  or  a significant  part  of  the  behavior  for  all 
the  data  displayed. 

The  derivation  of  the  magnetostatic  solution  is  explained 
in  the  next  section,  and  the  low  frequency  error  factor  is  the 
ratio  of  the  magnetostatic  solution  at  4>  = 0 for  a ^ b divided 
by  the  evaluation  for  a = b.  The  significance  of  this  low 
frequency  error  factor  is  also  discussed  in  the  next  section. 
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SFCTION  8 

APPROPRIATE  MAGNETOSTATIC  SOLUTIONS 

In  this  section  we  will  present  a solution  for  the 
magnetostatic  current  density  induced  on  an  ellipsoid.  By 
appropriately  choosing  the  parameters  describing  the  ellipsoid, 
we  will  obtain  the  magnetostatic  solutions  used  to  interpret 
the  finite  circular  cylinder  and  finite  elliptic  cylinder 
results.  The  parameters  of  the  ellipsoid,  the  orientation  of 
the  external  magnetic  field,  and  our  coordinate  system  are 
displayed  in  Figure  8.  We  choose  not  to  present  the  details 
of  our  derivation,  as  they  will  be  presented  shortly  in  a 
forthcoming  AFWL  technical  report  (7),  and  instead  only  present 
our  results.  The  current  density  at  points  in  the  z = 0 
plane  only  has  a z-component  which  is 

K(<t>,z  = 0)  = 2HQb  cos  <J)  (2  - aQ)  1 ( a 2 sin2  ^ + b2  cos2  <j> ) z 

where  x = a cos  <p 
y = b sin  <{> 

aQ  = 2a (a  + b) _1  (1  - 6) 

6 = b (a  - b)~1(l  - (b/c)2)"1/2(E(Y,0)  - (1  - (b/c)2)1/2) 

and  E (y , 8)  is  an  elliptic  integral  of  the  second  kind  requiring 
the  further  definitions 

E(y,B)  = J\  1 - 82  sin2  t)1/2  dx 

o 

where  02  = (1  - (a/c)2)/(l  - (b/c)2) 

Y = arccos  (b/c) 

When  a = b,  then  K(<f>,z  = 0)  reduces  to  the  spheriod  solution 
given  in  Reference  1.  In  Reference  1,  it  was  in  effect  stated 
that  a quantity  equivalent  to  our  definition  of  6 was  evaluated 
as  0.0203  for  a/c  = b/c  = 0.1.  This  ratio  of  0.1  was  chosen 
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to  correspond  to  the  length  to  diameter  ratio  of  the  circular 
cylinder  for  which  we  presented  our  data  in  Figure  6.  The 
value  of  <S  for  a/c  = 0.15  and  b/c  = 0.03  corresponding  to  the 
finite  elliptic  cylinder  data  presented  in  Figure  7 was  also 
investigated.  Due  to  the  fact  that  we  had  limited  tables  of 
elliptic  integrals  and  the  fact  that  our  intent  is  only  to 
justify  the  approximation  5 = 0 , we  are  content  to  say  it  is 
less  than  0.01. 

Finally  we  mention  that  6 does  approach  zero  as  c 
approaches  infinity.  The  formula  for  6 = 0 is 

2.2  2 2 -l/?~ 

K(4>,z  = 0)  = (a  + b)  H cos  <J>  (a  sin  $ + b cos  (p)  / z 

(1) 

and  this  solution  is  valid  whether  or  not  a and  b are  equal. 

For  completeness,  we  exhibit  the  case  where  a = b,  or 

K (4> , z=0 ) = 2H  cos  <}>  z (2) 

As  shown  in  Figures  6 and  7,  Equations  1 and  2 with  the  nor- 
malization H =1  describe  the  dominant  behavior  of  the  current 
o 

density  induced  in  the  vicinity  of  the  center  of  a finite 
elliptic  and  circular  cylinder  for  a band  of  frequencies. 

Having  established  this,  a significant  point  to  be  made 
by  this  paper  is  related  to  the  fact  that  Equations  1 and  2 
are  not  identical.  The  reason  we  make  this  statement  is  that 
if  long  wavelengths  justified  the  gross  modeling  of  physical 
features  of  structures,  then  one  would  not  expect  a significant 
difference  of  the  current  induced  on  a cylinder  having  an 
elliptic  cross-section  if  the  cross-section  were  modeled  as 
a circle.  A quick  way  to  obtain  a quantification  of  the 
error  is  to  compare  the  value  of  the  current  density  at  4>  = 0. 
The  resulting  current  density  for  a perfect  measurement  or 
calculation  on  a circular  cylinder  is  2 while  if  the  real 
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cylinder  had  an  elliptic  cross-section  the  value  should  be 
1 + a/b.  Relating  this  result  to  the  modeling  of  the  wing 
of  an  aircraft  where  we  might  choose  a/b  = 8,  then  we  see 
that  there  is  an  error  factor  of  9/2,  where  old  long  wavelength 
arguments  would  lead  one  to  believe  there  should  be  no  error 
at  all.  The  full  implications  of  this  result  go  far  beyond 
simple  finite  cylinders.  New  emphasis  in  deciding  the  amount 
of  detail  required  for  the  modeling  of  physical  features  of 
metallic  enclosures  is  warranted. 

Finally,  as  a secondary  issue  and  because  of  the  fact 
that  we  have  Equations  1 and  2 available  for  intrepretation, 
we  are  in  a position  to  clarify  a misconception  regarding  the 
magnetostatic  solution  for  complex  structures.  It  is  that 
2az  • (n  x H^)  happens  to  equal  the  magnetostatic  solution 
Equation  2;  however,  it  is  not  a universal  formula  as  can  be 
seen  from  the  fact  that  it  does  not  equal  Equation  1. 
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SECTION  9 


QUALITATIVE  LOW  FREQUENCY  CURRENT  DENSITY  DISTRIBUTION 
ON  A COMPLEX  AIRCRAFT  MODEL 

In  describing  this  low  frequency  distribution  we  refer  to 
Figure  9.  The  total  extent  to  which  we  want  to  describe  this 
distribution  is  contained  in  the  two  equations 

• A.  A 

t • K (r)  = -t'  • K (r ' ) (3) 

s • K (r ) = s'  • K (r ' ) (4) 

The  derivation  of  these  equations  will  appear  in  a forthcoming 
AFWL  technical  report  (7).  In  order  to  explain  these  equations.it 

AAA  A 

is  necessary  for  us  to  still  define  s,  t,  s',  and  t'  after 
noting  that  K(r)  and  K(r')  represent  the  surface  current  density 
evaluated  at  symmetric  points  with  respect  to  the  xy-plane. 

The  vectors  s and  s'  as  well  as  t and  t'  are  vectors  tangent 
to  the  surface  at  the  symmetric  points  r and  r ' . In  order  for 
the  relations  (3)  and  (4)  to  be  valid,  these  tangent  vectors 
must  have  a more  restrictive  definition  than  the  usual  one 

AAA 

which  is  that  (s,  t,  n)  form  an  orthonormal  set  on  the  surface 

A 

where  n is  the  outward  normal.  In  addition  to  this  requirement, 
we  will  describe  how  s is  defined  so  that  t will  then  be  deter- 
mined by  using  the  stated  definitions  of  n and  the  orthonor- 
mality to  give  t = n x s.  To  describe  s we  imagine  the  body 
cut  by  planes  parallel  to  the  yz-plane  and  define  s as  the 
tangent  vector  to  the  curve  formed  by  the  intersection  of  these 
planar  cuts  with  the  body.  In  Figure  9 we  display  the  curves 
and  s vectors  associated  with  two  cuts  through  the  body.  In 
this  figure  we  show  several  t vectors  that  are  convenient  to 

display  as  well  as  quantities  t and  s which  are  the  xz 
„ „ — xz  — xz 

projections  of  t and  s. 
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Now  that  the  meaning  of  f3)  and  (4)  has  been  explained, 
we  would  like  to  briefly  describe  the  utility  these  equations 
can  have  in  addition  to  their  providing  some  qualitative  insight 
concerning  the  induced  current  density  distribution.  We  first 
note  that  (3)  alone  is  sufficient  to  identify  the  error 
associated  with  current  density  calculations  obtained  by 
dividing  a bulk  current  by  a local  circumference.  This 
equation  also  provides  enough  information  to  draw  the 
conclusion  that  t • J<(r)  = 0 along  the  curve  formed  by  the 
intersection  of  the  xy-plane  with  the  body. 

Before  concluding  this  section,  we  will  present  a brief 
outline  concerning  the  derivation  of  (3)  and  (4).  This 
derivation  is  a generalization  of  the  derivation  presented  in 
Section  2 of  Reference  1.  It  is  possible  to  derive  an 
integral  equation  for  a quantity  that  when  set  equal  to  zero 
implies  (3)  and  (4).  It  is  argued  that  this  quantity  should 
be  set  equal  to  zero  in  the  low  frequency  limit  because  the 
source  term  for  the  integral  equation  vanishes  in  this  limit 
and  the  integral  equation  has  a unique  inverse.  Finally,  it 
can  be  shown  that  (3)  and  (4)  are  compatible  with  a magneto- 
static solution  for  the  aircraft  model  immersed  in  a static 
magnetic  field  oriented  in  the  y direction. 
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SECTION  10 
RECOMMENDATIONS 

It  has  been  demonstrated  that  for  an  extended  band  of  low 
frequencies,  the  induced  current  density  on  several  objects 
behaves  predominantly  like  a magnetostatic  solution.  We 
conclude  that  the  magnetostatic  response  of  an  object  should 
be  considered  in  determining  low  frequency  modeling  requirements 
for  an  external  coupling  analysis  or  measurement  program.  An 
appreciation  concerning  the  significance  of  magnetostatics  on 
determining  modeling  requirements  can  be  obtained  by  noting  the 
exact  analogy  between  magnetostatics  and  irrotational  and  incom- 
pressible fluid  flow  around  a rigid  body.  For  a rigid  perfectly 
conducting  body,  the  velocity  flow  lines  and  the  magnetic  field 
lines  in  the  vicinity  of  the  body  are  identical..  If  certain 
features  of  aircraft  or  missiles  would  have  a significant  effect 
on  fluid  flow,  then  they  would  have  a significant  effect  on  the 
current  induced  by  an  EMP . One  should  become  concerned  with 
modeling  such  features  as  missile  tips  and  fins  as  well  as  air- 
craft features  such  as  engines,  wing  cross  sections,  and  extended 
junctions,  particularly  when  the  point  of  entry  is  in  the 
proximity  of  these  features. 

Finally,  it  should  be  emphasized  that  magnetostatic  con- 
siderations should  only  be  part  of  a study  to  determine  modeling 
requirements.  The  primary  significance  of  magnetostatic  related 
requirements  is  that  tney  are  much  more  severe  than  what  has 
previously  been  thought  justifiable  as  a result  of  considering 
long  wavelength  far  zone  scattering  results. 

In  addition  to  the  recommendations  concerning  modeling 
requirements,  we  would  like  to  make  the  following  recommendations 

• A thorough  study  should  be  initiated  to  determine 

the  extent  to  which  error  and  confidence  requirements 
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for  external  coupling  analysis  should  be  diminished 
due  to  other  uncertainties  as  well  as  excess  hardening 
capabilities . 

In  the  interest  of  validation,  the  results  of  other 
measurement  programs  or  analyses  should  be  compared 
to  the  Harvard  group's  data  and  our  MFIE  code's 
calculations . 
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